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Abstract. In this paper, we analyze the conditions for convergence toward General 
Relativity of scalar-tensor gravity theories defined by an arbitrary coupling function 
a (in the Einstein frame). We show that, in general, the evolution of the scalar field 
(ip) is governed by two opposite mechanisms: an attraction mechanism which tends to 
drive scalar-tensor models toward Einstein's theory, and a repulsion mechanism which 
has the contrary effect. The attraction mechanism dominates the recent epochs of 
the universe evolution if, and only if, the scalar field and its derivative satisfy certain 
boundary conditions. Since these conditions for convergence toward general relativity 
depend on the particular scalar-tensor theory used to describe the universe evolution, 
the nucleosynthesis bounds on the present value of the coupling function, cto, strongly 
differ from some theories to others. For example, in theories defined by a oc| ip \ 
analytical estimates lead to very stringent nucleosynthesis bounds on a?o (<; 10 -19 ). 
By contrast, in scalar-tensor theories defined by a oc tp much larger limits on ao 
(< 10~ 7 ) are found. 
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1. Introduction 

Scalar-tensor (ST) gravity theories |L], 0, [| have become a focal point of interest in 
many areas of gravitational physics and cosmology. They provide the most natural 
generalizations of General Relativity (GR) by introducing an additional scalar field, 0, 
the dynamical importance of which is determined by an arbitrary coupling function 
Ld(4>). Indeed, most recent attempts at unified models of fundamental interactions, i.e., 
string theories |@ predict the existence of scalar partners to the tensor gravity of GR, 
and would have ST gravity as their low energy limit. In addition, ST theories are 
important in cosmology because they provide a natural (non-fine tuned) way of exiting 
the inflationary epoch. 

Solar system experiments (time delay in the Viking Mars data, Lunar Laser 
Ranging, etc.) can put limits on the present deviation of ST theories with respect 
to GR. It is admitted that such experiments impose the constraint ujq > 500 to the 
present value of the coupling function and, therefore, that ST theories are at present 
very close to General Relativity. This limit on u does not necessarily imply that the 
universe evolution is, at any time, very close to that found in GR. It has been shown 
H 0] that, in some ST theories, the cosmo logical evolution drives the scalar field toward 
a state indistinguishable from GR. Within this 'attraction mechanism', the scalar field 
can play an important role only in early cosmology because, afterwards, it evolves toward 
a state with a vanishingly small scalar contribution. 

The dynamical importance of the scalar field in the early universe can be checked 
by means of the observed abundance of light elements, which has to be explained as an 
outcome of the primordial nucleosynthesis process (PNP). Contrary to the weak field 
limit, there is not a commonly accepted PNP constraint on c^o- Some authors [§, [J 
have recently found that, in the framework of some ST theories, Uq ^ 10 7 is required 
to obtain the observed primordial abundance of light elements. Other authors [|l(], |l I 
have instead found that the PNP test typically imposes the limit uj ^ 10 20 . 

In this paper we will first reexamine the attraction mechanism of Refs. |6|, [jj and, 
then, we will investigate the reason for the enormous discrepancy (thirteenth orders of 
magnitude) in the PNP bound on ujq obtained by different authors. We will show that, 
in general, the evolution of the scalar field is governed by two opposite mechanisms: 
an attraction and a repulsion mechanism. The attraction mechanism dominates the 
recent epochs of the universe evolution only if the scalar field and its derivative satisfy 
certain boundary conditions which depend on the particular scalar-tensor theory used 
to describe the universe evolution. We will also apply this generalized formalism to 



the theories considered in Ref. ||11|| , where the coupling function was restricted to be a 
monotonic function of time. We will show that the nucleosynthesis bounds numerically 
obtained in |TT| (t^o ^ 10 20 ) are in close agreement with the analytical estimates for 
these theories. When the same arguments are applied to other ST theories (as those of 
Refs. H 0]), one obtains much less stringent bounds on ujq. Consequently, the particular 
ST theory used to describe the universe evolution has a crucial importance on the PNP 
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limit on ojq and, therefore, it is not possible to establish a general and unique limit for 
all ST models. 

The plan of the paper is as follows. We begin outlining the scalar-tensor theories 
as well as the two frames usually considered to build up cosmological models in their 
framework (Sec. 0). An autonomous evolution equation for the scalar field is then 
obtained in Sec. |3] both for the Jordan and the Einstein frame. Using this equation, 
and the particular family of scalar-tensor theories specified in Sec. f|, we then analyze 
the evolution of the scalar field both in the radiation-dominated epoch (Sect. f|) and 
the matter-dominated epoch (Sect. ||). Our results are then applied to estimate the 
nucleosynthesis bounds on ujq for this family of theories (Sec. [?]). Finally, conclusions 
and a summary of our results are given in Sec. || 



2. Scalar-Tensor Gravity Theories 

2.1. The Jordan frame 

We consider a scalar-tensor gravity theory in which the gravitational interaction is 
carried by the metric g^ u and an additional massless scalar field, <fi. Since these fields 
are measured through laboratory clocks and rods, which are made of atoms and are 
essentially based on a non-gravitational physics, the action describing a scalar-tensor 
gravity theory must keep unaltered the basic laws of non-gravitational interactions (as, 
e.g., statistical physics and electrodynamics). When this fact is taken into account, 
the resulting space-time units are termed as "atomic", "Jordan" or "physical" units 



T2L 13, 14, 15 



Using the Jordan frame, the most general action describing a massless scalar-tensor 
theory of gravitation is |], 0, || 

s = -i- fw- ^-<t>,»<f>>nV=gd 4 x + s M (i) 



16vr 

where 7Z is the curvature scalar of the metric g^ u , g = det^g^), is the scalar field, and 
u)(4>) is an arbitrary coupling function. 

The variation of Eq. (P with respect to g^ v and <fi leads to the field equations: 



g^O(f)) (2a) 



(3 + 2c;)n0 = 8nT - ^(f) a (f)' a (2b) 

dip ' 

which satisfy the usual conservation law 

= (3) 
where T^ u is the energy-momentum tensor and D0 = g^ v '<f> yfJi - u . 
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For a homogeneous and isotropic universe, the line-element has a Robertson- Walker 
form and the energy-momentum tensor corresponds to that of a perfect fluid. The field 
equations (Pal) and (EE) then become 



8tt _ c 2 K R 2 to $ 2 m 

3$ P ~ ~W + P 2 ~ e $2 + ^ ( a J 



(p + 3P/c) = 2 1 1 1 (46) 

3$ V ^ 11 R R<f> <f> 3 $ 2 

$ + 3-$ = - 1 [87rG(p - 3P/c 2 ) - ^$ 2 ] (4c) 
P (3 + 2cj) l VP 7 ' d$ J v 7 

where K = 0, ±1, $ = G0, P(i) is the scale factor, p and P are the energy-mass density 
and pressure, respectively, and dots mean time derivatives. In addition, we have the 
usual conservation equation: 

d{ P R 3 ) + {P/c 2 )dR 3 = (5) 

which ensures that the standard laws of non-gravitational interactions are not modified 
by the presence of a scalar field. 

2. 2. The Einstein frame 

When the metric is assumed to be measured through purely gravitational clocks and 
rods, the space-time units are termed as "Einstein" or "spin" units. In this frame, the 
general action describing a massless scalar-tensor theory can be obtained from Eq. (|T|) 
by a conformal transformation 

9i» = a2 {v)9% (6a) 
A 2 ( v ) = (6b) 

where A((p) is an arbitrary function related to co>($) by 

a = (3 + 2^/ 2 = ^ (7) 

dip 

Using Eqs. ([I]), (ffiiD, ( pty and (|7|), one obtains 

S * = I ( n * ~ 2W)V=9~*— + S* M (8) 

107rCr* J c 

where G* is Newton's constant and asterisks denote quantities expressed in Einstein 
units. Since our measures are based on non-purely gravitational rods and clocks, 
quantities written in the Einstein frame are not observable. Comparison between theory 
and observations must be then performed by using the Jordan frame. 
From the action (|8p, the Einstein field equations are: 

11% = 2<p,w iV + 8nG*(T; u - l -Tg;j (9a) 
□*p= - 4irG*a((p)T* (9b) 

— 1/2 

where T* = 2g* 5S m /5g% is the energy-momentum tensor in Einstein units. 
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If we consider a homogeneous and isotropic universe, the field equations (|9 
become 

- 3 ^ f S- toG «("« +3p « )+2 (S 2 (ioa) 

^hikwr-^'^'-^ (10c) 

and the 'conservation' equation is modified to 

d{p*Rl) + P*d(Rl) = (p, - ?>P*)Rlda{ip) (11) 

Since the mass-energy is not conserved in Einstein units, the basic laws of non- 
gravitational physics are modified in this frame (see |16], [17|, [18|] for a reformulation of 
nuclear reaction rates and thermodynamics in Einstein units). 

3. Decoupled evolution of the Jordan scalar field 

In the form given above, the time evolution of the scale factor and the scalar field are 
coupled both in the Jordan and the Einstein frame. Previous works || [7| have found 
that, by introducing an appropriate change of variables, it is possible to find an evolution 
equation for the Einstein scalar field which is independent of the cosmic scale factor. 
We will now show that it is also possible to find such a decoupled evolution equation 
for the Jordan scalar field. 

Let us define the functions 

V> = -ln$ (12a) 

7 = ^ (126) 
P 

2>c 2 K^ 



- S«G P R* (12C) 

W = (3 + 2u)/3 (12d) 

The Jordan evolution equations of R and $ then become: 

| + jfi + $ + 2W^ 2 = -^(1 + 37) (13a) 

^(l- e) =g + ^ 2 -^ (136) 

•■ R ■ 8%Gp . 1 dW ; , . . 

2 ^ + 4 ^ + 6-^=3^(l-3 7 )--— r (13c) 

In order to obtain a decoupled evolution equation for if), we will define a 'time' 
parameter, p, as: 

R 

dp = h c dt; h c = - + tp (14) 
ti 
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In terms of these variables, and denoting /' = df/dp, Eqs. ( |13«| )-( |l3c|) reduce to: 

^^-l-2^-^(l + 3 7 ) (15) 



^(l-^' 2 ) = ^(l- e ) (16) 

ip" + - iP' 2 + Zip' = 
h c 

47rGpl-3 7 1 dW 2 

and, using Eqs. (|i~6|) and ( |l5l) to eliminate /i c and /i' c , respectively, we finally obtain: 
2(1 -e) 



-ip" + (3 - 3 7 - 4e)</>' 



1 - Wip' 2 

l-3 7 / 1 - e \ 1 ^ „, 

-ip 



W \ 1 - W^' 2 / W dip 

which is an evolution equation for ip = (1/2) In $ independent of the evolution of the 
cosmic scale factor. This equation is the analogous, in the Jordan frame, to Eq. (3.15) 
of ||. Therefore, we find that a decoupled evolution of the scalar field is not an exclusive 
feature of the Einstein frame. 

In order to find the well-known decoupled evolution equation for the Einstein scalar 
field, we will first multiply Eq. (H) by W^ 2 (l - Wip' 2 )/ '2(1 - e). We obtain: 



Wip 



12 



' J>7 -3(1-7- U)W 1/2 i>' (19) 



W 1 ' 2 v ' 3 

We can now change to the Einstein frame by performing the conformal 
transformation Q6fl|)-(|6l]) with 

^ = _^y = _ a -i (20) 
dip 

Eq. ( pj|) then becomes: 

^—4<P" + (1 - 7 - Uv' = -«(1 - 3 7 ) (21) 
3 — ip' z 3 

which agrees with Eq. (3.15) of 0. 

We also note that, by defining if* and p* through 

= hc . p = & pmj (22) 

Eq. (HD yields 

SvrGp, _ 3-y/ 2 

" t^t (23) 

From this equation we see that, when e = 0, the local positivity of the energy 
density implies that 

y?' 2 < 3. (24) 
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Most papers analyzing the convergence toward General Relativity of scalar-tensor 
theories are based on the Einstein frame. In order to make easier the comparison of our 
results with those founds in previous works, we will use hereafter the Einstein frame. 
Nevertheless, it must be noted that most of our conclusions can also be found by using 
the Jordan evolution equation (|T8|). 



4. The coupling function 



The time evolution of scalar-tensor theories can be studied only after specifying a 
functional form of the coupling function. Barrow and Parsons |19| have noted that 
most expressions used in the literature for (3 + 2u) can be classified into three different 
families of theories: 

1 



Theories- 1: 3 + 2uo 



b) Theories-2: 3 + 2uj 



c) Theories-3: 3 + 2uj 



1 



Silln^l 5 
1 



II 



where Bi is an arbitrary positive constant. 

The three theories defined by Eqs. (|25q|) -(|25 



(5 > 1/2) 
(5 > 1/2) 
(5>0) 



(25a) 
(256) 
(25 c) 



the early universe, which have been analyzed in detail by Barrow and Parsons [19 



imply very different behaviours of 

The 



first class of theories has also been studied by Garci'a-Bellido and Quiros ||20|| , Serna and 
Comer et al. 



Alimi [21 



22 



23 1 . However, it is important to note 



and Navarro et al. 

that all these theories have similar behaviours in the limit close to General Relativity 
($ — > 1). As a matter of fact, since ln$ ~ $ — 1 and $' 5 — 1 ~ 5 In $, we can take 

1 1 



3 + 2u 



(5 > 1/2) 



(26) 



£i|ln$| 5 B^i/jI 5 

to represent the way in which these three types of theories approach the limit of GR. 
Using the coupling function given by Eq. fl26|) , integration of Eq. (|20| ) yields 

sign(0) 



5)B l/2 



|2Vl 



(2-<5)/2 



(5<2) 



(2 - u,u l 

where we have normalized the integration constant so that if 



(27) 



corresponds to $ 



(or ■0 = 0). Note that, according to Eq. 

By introducing Eq. (|27|) into Eq. (| 
function: 



where B 2 



a = ^Iv 9 ! 2 " 5 — K (^)M 
£i/M( 2 _ 6 y/(2-s) >0) and 

2(6-1) 

K(<p) = B 2 \ip\ a-« 



3)» si & n (f) = -sign(0). 
we obtain the Einstein form of the coupling 



(2f 



(29) 
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5. Radiation-dominated evolution 

In the radiation-dominated epoch, the state equation is P/c 2 — p/3 (i.e., 7 = 1/3) and 
the curvature effects are negligible (e = 0). Consequently, the evolution equation of the 
Einstein scalar field (Eq. ^1|) is well approximated by: 

if" + -if' = (30) 



3-f' 2r 3 

which does not depend on the functional form of ct((p). 
The integration of Eq. (|30|) gives: 



where k is related to the initial (p = 0) velocity <f' R through: 

k 2 = ^fry 2 (32) 
3 - Wr) 2 



In terms of k, the solution of Eq. (^) is: 
V 9 = <Pr : ~ sign (A;) In 



VI + k 2 e- 2 P + ke~ p 



VTTW+k 

6. Matter-dominated evolution: The attraction-repulsion mechanism 



(33) 



Let us now analyze the evolution of the scalar field during the matter-dominated era 
(7 = 0) of a flat universe (e = 0). The evolution equation for the Einstein scalar field 
is, in this case: 

2 + (p' + a = (34) 



3- cp' 2 

As in previous works || |§, we will first assume that, at some time (for instance, 
at the beginning of the matter-dominated era), the scalar-tensor theory is not very far 
from GR so that we can neglect if' 2 against 3 and, in addition, the coupling function 
a{ip) is represented by Eq. (f2~8|). 

In this case, the evolution equation (133) reduces to: 



where 



2 

-<f" + f' + <r v K((p)(p = (35) 



a v = sign(y?) (36) 

When = +1, the above expression corresponds to the evolution equation of 
a damped harmonic oscillator with a variable elastic coefficient n(f). The first term 
(2(f"/3) represents the total force on a fictitious particle of mass m = 2/3. The second 
term (f') corresponds to a friction force proportional to the velocity and, finally, the 
third term represents an elastic force with a variable coefficient n(f). The existence, 
in these conditions, of a damped oscillatory behaviour of the Einstein scalar field was 
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first reported by Damour and Nordvedt M and it is usually termed as the 'attraction 
mechanism' toward General Relativity . 

We note however that, if a v = —1, the effective elastic coefficient a^n^tp) is negative 
and, consequently, there exists a 'repulsion' mechanism instead an attraction one. We 
will now analyze the matter-dominated evolution of the scalar field by considering 
separately the cases 5 = 1 and 5 ^ 1 in Eq. fl2"9"|). 

It must be noted that the class of scalar-tensor theories analyzed in || correspond 
to a positive constant a v and, therefore, they are always attractive. 

6.1. Case of a constant elastic coefficient (5 = 1) 

When 5 = 1, the elastic coefficient defined by Eq. (^) becomes a positive constant B 2 . 
The scalar field evolution equation then reduces to: 



the general solution of Eq. fl3"T|) admits four different behaviours, depending on B 2 and 



<p" + -<p' + -cr^if = 



(37) 



Since the roots of the characteristic equation are: 




(38) 



0", 



a) Damped harmonic motion (a v > 0, B 2 > 3/8) 




(39) 



where 




(40 a) 



(406) 



(40 c) 



b) Critically damped motion (<r v > 0, B 2 = 3/8): 




(41) 



where 



Ci = (f'o + -v?o; C 2 = v? 
c) Overdamped attraction motion (a v > 0, B 2 < 3/8): 



(42) 



(43) 



where 
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A = -yi--S a <- (44a) 
Ci = ^0 + ^ + 3/4)^. 



2/5 



C 2 = - ^o + (-^ + 3/4)^0 (44c) 
d) Overdamped repulsion motion (a v < 0, and any B 2 ): 



3 



<p = e - ! P[C 1 eP* p + C 2 e- p2P ) (45) 



where 



Ci = ^o + (3/4 + ^ 0; 




2/5 ; 



2 



C 2 = _ rt + (3/4-A )Wl 

Apparently, only the first three solutions converge toward GR (<p = 0) while the last 
one diverges to ±oo as p — > oo. However, it must be noted that each of such solutions 
governs the whole matter-dominated era if, and only if, a 9 does not change. This is 
the case of the models analyzed in Refs. || where, in Eq. (|37f), is forced to be 
a v = +1 at any time. 

In the case of the models considered here, where the sign a v of <p can be variable in 
time, the matter-dominated evolution of the scalar field has a much more complicated 
behaviour. As a matter of fact, the above four solutions (Eqs. |39|-[45|) admit the 
possibility of a change in the sign of ip at a finite time p > (this is specially obvious for 
the first solution due to the oscillatory behaviour of the cosinus function). Therefore, 
an initially convergent model (a^ = +1) could finally diverge from GR as a consequence 
of a change in the sign of ip. Reciprocally, an initially divergent model (a^ = —1) could 
finally converge toward GR if sign(<y?) changes. We must then perform a more detailed 
analysis to find the conditions for convergence (or divergence) with respect to GR. 

According to Eqs. (|39|)- (|46c|) , ip vanishes (and, therefore, a v changes) at a finite 
p > if, and only if, the initial values of ip and ip' satisfy the following conditions: 

a) always 



b ) <Po < -f^o, (<A) > 0) 

c) <f'o < -(f+/3i)<A), (<A)>0) 



(47) 



d) <p' >-(l+^)vo, (<A)<0) 

When these conditions are satisfied, an initially attraction behaviour becomes a 
repulsion one at p\ > 0. Reciprocally, an initially repulsion behaviour becomes, at 
Pi > 0, an attraction one. The question is now to determine whether these new 
behaviours, reached at p > pi, govern the rest of the matter-dominated evolution. 
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Let us first consider the attraction-to-repulsion case. In this situation, the 
attraction models given by Eqs. fl39]), ( pLTj ) and ( pE5|) (with p < p x ) imply a vanishing 
scalar field at p = p±. After this time, the sign of ip has changed and, therefore, all 
models become expressed by Eq. (f4"5|) but with the integration constants (see Eqs. [46 1\ 



and [46 dj given by: 



Ci = -C 2 = ^- (48) 
where denotes the ip' value at p = pi, and we have taken into account that <p(pi) = 0. 



Using Eqs. (45) and (|4~8j) , the scalar field evolution for p > p\ is given by: 



Ci 


_ WW 


c 2 






Ci 


= tfl, 


c 2 


Ci 


_ vi 

2ft ' 


c 2 



cp = I±. e -\v[eP>* - e -^} (49) 
where p = p — p\. 

This solution never vanishes for p > p\ and diverges to ±oo when p — > +oo. 
Consequently, none of these models finally converges towards GR. 

Let us now consider the repulsion-to-attraction case. In this situation, the repulsion 
model given by Eq. (f4~5|) (with p < pi) implies a vanishing scalar field at p = pi. 
After this time, the sign of ip changes and, therefore, the scalar field evolution becomes 
expressed by Eqs. fl3~9|), ( |4l]) or (|43|) , depending on the B 2 value. The integration 
constants in this new regime are: 

tt/2 [B 2 > 3/8) 

(B 2 = 3/8) (50) 
-Ci {B 2 <3/8) 

Therefore, when p > p\. 

^e^fPsin(cJip) (B 2 > 3/8) 

<p{e-^p (B 2 = 3/8) (51) 

^ e -|p[ e ftP _ e -ftP] (b 2 < 3/8) 

We then find that, if B 2 > 3/8, the new attraction behaviour is oscillatory. 
Consequently, it will return later to have a repulsion motion similar to that of Eq. 
(|49|), and these models will finally diverge from GR. On the contrary, if B 2 < 3/8, Eqs. 
(pl|) imply that ip never vanishes at p > pi except for p — > +oo and, therefore, these 
models are finally attracted toward GR. 

Summarizing, all models with B 2 > 3/8 diverge from GR, whatever the initial 
values of a^, ip and ip' are. In the same way, models with B 2 < 3/8 and a v > (at 
p = 0) will diverge from GR if conditions (f4~7|)b or (f47|)c are satisfied, while models 
with a v < (at p = 0) will diverge from GR if conditions (|7)d are not satisfied. 



Consequently (see Fig. 1), the scalar-tensor theories defined by Eq. (|28|) with 5=1 will 
converge toward GR if, and only if B 2 < 3/8 and: 

^ ~ (3/4 + Pi)<po (for any a v value) (52) 
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Figure 1. Matter-dominated evolution of <p in converging models with tpo < 
and ip' = ificrit (solid line), ipa < and tp' > (p cr it (dotted line), (po > and 
tp' Q = tpcrit (dashed line), and (fo > and <p' a > (^ cr it (dotted-dashed line), where 

<Pcrit = -(3/4 + ^)^0- 



This result can be understood by considering the scalar field (p as the position 
of a particlelike dynamical variable which moves under a potential V((p). Since the 
right-hand side of Eq. ([H]) represents a force term, this potential is given by: 

V(<p) = (1 - 3 7 ) J a{y)dy (53) 

If oe(ip) is given by Eq. (|28| ) with 8=1, the matter-dominated form of V(ip) 
becomes: 

= \cr^ (54) 
where we have taken the origin of potentials so that V = at (f = 0. 
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-4-2 2 

Einstein scalar field 



Figure 2. Potential V(ip) for scalar-tensor theories with 6=1. Here, the particular 
case Bi = 1/9 is shown 

Figure 2 shows this potential. We see that V(<p) has a stationary point at <p = 
which does not correspond to a local minimum or maximum. Particles with <p > are 
attracted towards cp = while those with cp < are rejected from <£> = 0. 

In the case of ST theories characterized by an inefficient friction [B 2 > 3/8), any 
particle will 'slide' on this potential and will move toward ip — > —00, whatever the initial 
values of <p and ip' are. This explains our previous result that any theory with B 2 > 3/8 
diverges from GR. 

In the opposite case of theories characterized by an efficient enough friction (B 2 < 
3/8), particles with ip > have an overdamped or critically damped motion. Therefore, 
if the initial velocity ip' is not more negative than a critical value — (f3\ + 3/4)</? > the 
energy is completely dissipated by the friction term, and the particle stops at cp — 0. On 
the contrary, if ip < 0, the tendency to diverge from GR (<p = 0) can only be overcome if 
the particle has an initial velocity high enough to reach a point placed in the attraction 
(ip > 0) region. This explains our previous result that, in ST theories with B 2 < 3/8, 
the convergence towards GR impose on ip' the constraints given by Eqs. ( [52]) . 

6.2. Case of a variable elastic coefficient (5^1) 

When 5^1, the matter-dominated potential is given by 

Vlp) = 2 -^a,B 2 \ V \^ (55) 

Since the form of this potential is qualitatively similar to that of Eq. (|54]), we can 
expect for these models the same types of behaviours as those found in the 5=1 case. 
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Nevertheless, as we have seen in the previous subsection, these behaviours strongly 
depend on the relative importance of the friction term and on the initial value of ip 
and ip' . We can obtain some insight on the conditions for convergence toward GR by 
considering the two following limiting cases. 

When 1 < 5 < 2, the elastic coefficient defined by Eq. ( p9|) tends to zero as ip — > 0. 
Consequently, in the limit close enough to GR, we can neglect the term a ip n{ip)ip of 
Eq.(|35|), so that the scalar field evolution equation becomes friction dominated: 

v " + ^' = (56) 

The solution of this equation is: 

V = Ci + C 2 e^ p = po + ^o(l - e~l p ) (57) 

which implies that ip vanishes at a finite time p% > provided that: 

3 

3 

< -7><Po (<^o > 0) (586) 

These expressions correspond to the limit B 2 — > (very strong friction) of Eq. (|4~7|) 
for the 5=1 case. Therefore, according to Eqs. fl52|), theories with 5 > 1 will converge 
toward GR provided that 

> ~l<Po (59) 

On the contrary, when 1/2 < 5 < 1, the elastic coefficient becomes infinity as <p — > 
(see Eq. ^9|). The friction term of Eq. ([H]) is then very inefficient and, therefore, ST 
theories behave as the limit B 2 — > c>o (very weak friction) of the 5 = 1 case. That is, 
the particle will 'slide' on the potential and will move towards ±oo. These models 
will then diverge from GR whatever the initial values of ip and ip' are. 

6. 3. Influence of the mass term 

In the preceding discussion, we have obtained necessary and sufficient conditions for 
convergence toward GR in the limit of small ip' velocities. We will now show that, in 
the most general situation with a non- negligible ip', there also exists a necessary (but 
not sufficient) condition for convergence toward GR. 

Let us consider Eq. fl34|) with the coupling function defined in Eq. (|28|): 
2 

-ip" + ip' = -a v K{ip)ip (60) 

3 — ip! 1 

By integrating this equation, one obtains 

inf#±^^r+(^ )= 



- b 2 r \ip\^dp (6i) 

J 
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Let us now assume that theories converge toward GR as p — > +00. In this case, 
<f — > and ip' — > as p —* +00, so that Eq. ([H]) becomes 



1 

5 



lnf V ^ + ^ ^+y = ^2 / M^dp>0 (02) 



v V3-^' y j o 

Therefore, the assumption of convergence toward GR lead to the necessary 
condition: 

<*> > ^TT^ (63) 

which leads to tp' > —3ipo/2 when ipo <C 1 while, in the limits tpo — > ±00, implies 
¥>o > Tl 

7. Implications on Primordial Nucleosynthesis 

The primordial nucleosynthesis process (PNP) starts, in the early universe, soon after 
the cosmic temperature becomes lower than that needed to maintain the proton-to- 
neutron ratio in its equilibrium value (freezing-out temperature). The nuclear reactions 
that then take place lead first to the Deuterium and Helium-3 formation. These elements 
are then burnt to produce Helium-4 and very small amounts of heavier elements. This 
process is sensitive to the universe expansion rate during nucleosynthesis, which can 
be parametrized by the speed-up factor £ = H/H GR , where H is the physical Hubble 
parameter and Hqr is that predicted by GR at the same temperature. Since scalar- 
tensor theories predict a universe expansion rate which differs from that obtained in the 
framework of GR, the light-element production can be used to place constraints on the 
£ value at the beginning of the PNP process, £pnp- 

Most of works agree on the PNP constraint on £pnp (subscripts PNP denote values 
at the beginning of nucleosynthesis). When a conservative choice for the observed 
primordial abundances is considered, (,pnp is constrained to be 0.8 ^ (,pnp ^ 1-2 
f9|, [24j| . However, when a more severe choice for the observed abundances is taken, the 



PNP limits on £,pnp are 0.95 ^ £pnp ^ 1-03 [pj , |25| . The problem arises when this 
bound is used to place constraints on the present value of the coupling function, u . 
Serna & Alimi [Tl| have considered a family of scalar-tensor theories where the 



matter-dominated potential V(ip) is similar to that of Fig. 2 and, in addition, the 
coupling function 3 + 2uo is a monotonic function of time (i.e., 3 + 2uj diverges to 



20 



infinity only in the limit of very large times). Except for some particular cases |26 
the nucleosynthesis bounds obtained on these theories are very stringent: u ^ 10 
Other authors @, || have instead considered a different family of theories where V(y) 
has a constant sign and 3 + 2a; is not necessarily a monotonic function of time. The PNP 
bounds on these last theories are about thirteenth orders of magnitude lower (uo ^ 10 7 ) 
than those obtained in the previous case. 

It has been suggested in the literature that the above discrepancy in the PNP 
bounds on uo could be perhaps due to the existence in Ref. |TTJ of numerical instabilities 
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(runaway solutions) arising when the field equations are solved from a backward time- 
integration. We will now show that the integration sense has no effects on the PNP 
bounds (see the Appendix A for a test on the numerical stability against the time- 
integration sense). The important point to explain the above discrepancies is the 
particular scalar-tensor theory used to describe the universe evolution. 

We will consider the same family of scalar-tensor theories as in Ref. and we 
will restrict our discussion to singular models with 3 + 2u; > and a monotonic evolution 
of the speed-up factor. 

Let us first consider the case $ > 1 (ip < 0), where < and the scalar field 
enters the matter-dominated era with the overdamped repulsion motion of Eq. fl45l) . 
According to Eq. fl52|) , the convergence toward GR requires an initial scalar field velocity, 
ip' M (subscripts M denote values at the beginning of the matter-dominated era), with a 
positive value: 



¥m (64) 



Therefore, since pu < 0, the choice (p' M = is forbidden for this family of scalar- 
tensor theories. In addition, the requirement of a monotonic time evolution for the 
coupling function eliminates the possibility of (p' M > — (3/4 + /? 2 )v 9 m, where convergence 
toward GR is obtained after a rather complicate behaviour of 3 + 2a; in which it becomes 
infinity (when p = for the first time), then it decreases until reaching a local minimum 
value (maximum value of ip), and finally it increases monotonically to infinity (see Fig. 
2). 

Taking the only possible choice for these theories 



v'm 



Vm, (65) 
Eq. ( pEo| ) reduces to: 

tp = ip M e-(i +f32)p (66) 

where, from Eqs. (|46fcQ - (|46(fJ , fi% has a constant value larger than 3/4 (very overdamped 
motion, B 2 0) and smaller than 3\/2/4 (critically damped motion, B2 — > 3/8). 

Using Eqs. (p8|) and (|66|) , the present value of the coupling function is (with 5 = 1): 

al = BW M e~^ + ^ (67) 

where po ~ 10 (see Ref. @) is the amount of p time elapsed since the beginning of the 
matter-dominated era. Taking the less stringent value B2 ~ 1/9, Eq. flBTD implies 

al « KTVm (68) 

while using the largest value B 2 = 3/8 

a 2 « lO" 17 ^ (69) 

In order to estimate al, we must specify a <pm value compatible with the observed 
abundance of light elements. To that end we note that, according to Eqs. (0) and 
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), the physical Hubble parameter H is related to that measured in Einstein units 



by 



H 2 



Hi 

A 2 



[1 + ot(p 



l\2 



Therefore, using Eq. fl23|) , we obtain 
9 8nGpA 2 , 
3 — y>' 2 
and the speed-up factor £ 
3A 2 



A2 



H/ Hgr (where #£ R 



8nGp/3) is given by: 



-,'2 



1 + cap ) 



/\2 



(70) 



(71) 



(72) 



3 — y> 

which, in the limit of very small (ft values, implies £ ~ A 

In models with a monotonic evolution of the speed-up factor, the compatibility 
with the observed abundance of light elements imposes 0.95 <^ £ <^ 1.03 just before the 
big-bang nucleosynthesis process. Since £ ~ A and if) = —InA, this constraint implies 
(see Eq. gTJ) : 

/ -0.72 (if B 2 = 1/9) 
^ PPN W \ -0.39 (if 5 2 = 3/8) 

Taking into account the condition (ESI) for convergence toward GR (i.e., (p' M 
at the end of the radiation era), Eq. ( |3*ID implies k rs — (fM^ PR / v3 and Eq. (| 
yields 



(73) 
then 



v^ln 



<^| f e 2 P« + 3 - <^ M e 



PA 



+ 3 - y? M 



(74) 



where p# denotes the p-time elapsed from the relevant epochs prior to the big bang 
nucleosynthesis up to the end of the radiation era. Using fppN values of Eq. ([73]) and 
the rather conservative value Pr = 3 (i.e., the cosmic temperature decreases in about 
three orders of magnitude) , Eq. (|74f) implies: 

-3.7 • lO- 2 (if B 2 = 1/9) 
10~ 2 (if B 2 = 3/8) 



-2.0 



(75) 



Consequently, 




(if B 2 
(if B 2 



1/9) 
3/8) 



(76) 



in good accord with the numerical computations of Serna and Alimi JTT 



The case $ < 1 ((p > 0) requires a different analysis. Since the models 3 + 2u > 
and $ < 1 studied in |TTJ are all nonsingular (class-3), the cosmic temperature 



has a maximum value, T max , which corresponds to the minimum of the scale factor. 
Consequently, in this class of models, the bounds on a are mainly imposed by the 
condition that T max must be high enough to allow for the existence itself of primordial 
nucleosynthesis. The numerical computations of [I I] lead to limits on a 2 , close to those 
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given by Eq. (|76|). This nonsingular behaviour is found for any ip' M value, except for 
the (fine tuned) choice (p' M = where the scalar field evolution is frozen to a constant 
value during the radiation epoch. 

It is worth commenting the fact that the strong limit expressed in Eq. ([75) refers 
to scalar-tensor models with a monotonic evolution of £ (T) . As we remarked in previous 
works ||11|| , these limits lead to cosmological models which do not significantly differ from 
the standard GR one except perhaps for very early epochs before nucleosynthesis. On 
the contrary, when the speed-up factor has not a monotonic behaviour during the big- 
bang nucleosynthesis, the compatibility with the observed abundance of light elements 
is possible even for very large £ppn values [ll], |26| . In the framework of this last family 



of models, the nucleosynthesis constraints on (fM and a 2 , are not very stringent and, in 
addition, the allowed range for the baryon density is much wider than in the standard 
GR cosmologies. 

8. Conclusions 

In this paper we have analyzed the convergence of scalar-tensor theories (ST) toward 
GR and its consequences on the nucleosynthesis bounds on the present value of the 
coupling function. 

To that end, we have deduced an autonomous evolution equation for the Jordan 
scalar field. By writing this equation in Einstein units, we have analyzed the evolution 
of the scalar field both in the radiation-dominated epoch and in the matter-dominated 



epoch. We have considered a coupling function defined by Eq. (|28|), which reproduces 



all the models proposed by Barrow & Parsons |B1 in the limit close to GR. We have 
then shown that, in general, the evolution of the scalar field is governed by two opposite 
mechanisms: an attraction and a repulsion mechanism. The attraction mechanism 
dominates the recent epochs of the universe evolution if, and only if, the scalar field 
and its derivative satisfy certain boundary conditions which depend on each particular 
scalar-tensor theory. 

Our results have been then applied to obtain an analytical estimate of the Big-Bang 
nucleosynthesis (BBN) bounds on u . We have found that the particular ST theory used 
to describe the universe evolution has a crucial importance on the BBN limits on ujq. 
Therefore, it is not possible to establish a general and unique limit for all ST models. In 
the particular case of the theories analyzed in this paper, where a oc \ip\, our analytical 
estimates are in close agreement with the nucleosynthesis bounds numerically obtained 
in [11] (q!q ^ 10~ 20 ). Theories different from those analyzed in this paper could imply 
very different BBN bounds. For instance, in the case of a ST theory defined by a(f) oc ip, 
where only the attraction mechanism is present, the BBN bounds are (a 2 , <^ 1CT 7 ) M, |J. 
In the same way, in scalar-tensor theories with a non-monotonic evolution of the speed- 
up factor, the BNB limits are drastically relaxed to a value comparable to that obtained 
from solar system experiments (a 2 , <; 0.02). In addition, in this last case, the allowed 



range for the baryon density is much wider than in the standard GR cosmologies flTl], [26[ . 
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-log„(T) 



Figure 3. Numerical solutions for £ (solid line), log 10 [l + a -2 ] (dotted line) and <p' 
(dashed line) in an open (k = — 1) scalar-tensor theory defined by a 2 — (A 2 /3)(l — $) 
with A 2 = 0.9. The results obtained from a backward and a forward time-integration 
are represented by lines and crosses, respectively. 
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Appendix A: Numerical stability of the time-integration sense 



We will here comment a series of computations which were performed in |TT|] to test 
the numerical stability against the time-integration sense. In these computations, the 
ST cosmological equations were solved in the two senses (forward and backward in time) 
by using a sixth-order Runge-Kutta method which takes the temperature as variable. 
A standard particle content (baryons, electron-positrons, photons and neutrinos) of the 
universe is assumed, so that transition phases are solved without approximations. 

Figure |3] shows a typical example of the numerical results found for the family 
of scalar-tensor theories considered in |2T], |TT[. The lines in this figure represent the 



solutions obtained for £ (solid line), log 10 (l 



a 



-2> 



(dotted line) and </?' (dashed line) 



from a backward time-integration. The final results of this integration were taken as 
boundary conditions to a reverse (forward in time) integration. The results obtained in 
this last case are represented by crosses in the same figure. 

As we can see from Fig. 3, both integration senses give almost exactly the same 
results. This fact can be very well understood if one analyses the propagation of 
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numerical errors during the integration. To that end, we consider the simple case of 
a critically damped motion (B 2 = §), where the time evolution of (p is given by Eq. 

(0): 



<p(p) = Vo 




(77) 



Let now suppose that we integrate this equation backward in time (p varying from 
today up to large negative values in the past) with numerical errors Aip <C <p and 
A(p' Q <C <p'q on the initial values. Obviously, since the functions (p and (p' increase in the 
past, their absolute errors will also increase. However the relevant quantity to estimate 
the accuracy of the integration is not the absolute error but the relative error which, at 
any p-time, is given by 



(p) 




A^o + (jAyp + A<p' )p 



In the limit where p goes to 



-oo, the above expressions become 




|Av9o + A^q 
1^0 + (f/ 



(78) 
(79) 

(80) 



Since A<p -C y^o and A<p' <C <p' , we then deduce that the relative errors remain very 
small in the past. This explains why the computation presented in Fig. 3 are stable in 
both integration senses. 
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